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The specific-heat jump AC at a critical temperature T c in an anisotropic superconductor contain- 
ing both potential and spin-flip scatterers is calculated within a weak-coupling mean-field approxi- 
mation. It is shown that the presence of even a small amount of spin-flip scatterers in the sample 
leads to a drastic change in the dependence of AC on T c in a disordered (d + s)-wave or a strongly 
anisotropic s-wave superconductor. The implications for experimental tests for the presence of an 
s-wave admixture in the superconducting order parameter of high-T c superconductors are discussed. 



I. INTRODUCTION 

Although there is a considerable evidence in favour of predominantly d-wave in-plane symmetry of the super- 
conducting order parameter A(p) in high-temperature superconductors (HTSCs) [Q, a number of experiments 
point to an admixture of the s-wave component to d-wave A(p), implying a mixture of d-wave and s-wave 
components, see references in Ref. pp, or even to a strongly anisotropic s-wave A(p), see Refs. The non- 

pure d-wave symmetry of A(p) is supported in part by a long tail suppression of the critical temperature, T c , 
by defects and impurities jUI^] characteristic for a nonzero value of the Fermi surface (FS) average (A(p))ps, 
while (A(p))ps — for a d-wave pairing. However, neither the presence nor the absence of an isotropic s-wave 
component of A(p) in HTSCs has not been demonstrated unequivocally yet. 

Recently Harah et al. 6 have shown that significant features attributed to the s-wave part of A(p) may 
be seen in the measurements of the specific heat jump AC at T c in disordered (d + s)-wave superconductors 
with nonmagnetic impurities. It was found in Ref. that there is a qualitative difference in the evolution of 
AC under disordering in d-wave and (d + s)-wave superconductors. While in a d-wave superconductor, the 
value of AC normalized by the normal state specific heat C n (T c ) monotonously goes to zero as T c is suppressed 
by nonmagnetic disorder, in a (d + s)-wave superconductor there is a minimum at the curve of AC/C„(T C ) 
versus T c . The location of this minimum depends on a specific weight of an isotropic s-wave component in the 
(d + s)-wave A(p), i. e., on a specific value of (A(p))^s- It was suggested in Ref. that this effect may be 
used as a test for the presence of an s-wave admixture in HTSCs. 

The authors of Ref. jfi] restricted their consideration to the case of nonmagnetic disorder (i. e., potential 
scattering of charge carriers) only. Note, however, that a lot of experiments give evidence for the presence of 
spin-flip scatterers (along with potential ones) in non-stoichiometric HTSCs, e. g., in oxygen-deficient, doped or 
irradiated samples (see, e. g., the references in Ref. [Jj). The spin-flip scattering of charge carriers by magnetic 
impurities (even though they may be present in small proportion) can strongly modify the value of T c as well 
as other characteristics of a superconductor in the vicinity of T c El - It is therefore instructive to elucidate if 
an account for spin-flip scattering will change the conclusion drawn in Ref. HJ concerning the impurity effect 
on AC in anisotropic superconductors. 

The influence of the ratio of spin-flip to non spin-flip scattering rates on the specific heat jump in a weakly 
anisotropic superconductor has been considered earlier by Okabe and Nagi P]. In the present paper we study 
the combined effect of potential and spin-flip scatterers on AC in a superconductor with arbitrary anisotropy 
of A(p), including s-wave, d-wave and (d + s)-wave as particular cases. In what follows, we make use of the 
following approximations: (i) we consider the pairing interactions in the weak-coupling limit of the BCS model 
and assume that the pairing potential can be taken in a factorizable form; (ii) we treat the impurity scattering in 
the Born limit; (iii) we assume an s-wave scattering of charge carriers by both potential and spin-flip scatterers, 
so that the scattering matrix elements are momentum independent; (iv) we neglect any dynamic pair breaking 
effects pJi|. 

One should keep in mind that the mean-field BCS theory does not describe the effects of spatial variation 
of the order parameter in the vicinity of impurities [H]. These effects may be significant in superconductors 
with short coherence length. However, experiments on the impurity induced suppression of T c in anisotropic 
superconductors are, in fact, rather well described within the mean-field Abrikosov-Gor'kov approach []2J (see, 
e - g-> So, the spatial variation of the order parameter may appear to have a little influence on the physical 

characteristics, including AC, in the vicinity of T c . Note also that strong-coupling corrections result just 
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in renormalization of the scattering rates ^Sj and overall increase in the magnitude of AC, without qualitative 
changes in the dependence of AC versus T c |16M17| . As for the use of the weak scattering (Born) limit, it seems 
to be justified by a very close similarity of AC versus T c curves calculated in the Born and unitary limits for the 
case of purely nonmagnetic disorder in d-wave and (e? + s)-wave superconductors [HIElEj- We set fi = Hb = 1 
throughout the paper. 



II. FORMALISM 

The specific heat jump is defined as AC = C S (T C ) — C„(T C ), where the subscripts s and n refer to the 
superconducting and normal state respectively, C„(T C ) = (2n 2 /3)N(0)T C , and N(0) is the density of electron 
states per spin at the Fermi level. The value of AC can be expressed in terms of the thermodynamic potential 
Q, as 



AC = -T, 

T=T, 

For the value of the difference fl s — tt n one has (see, e. g., Ref. \12\ ) 
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where Vq is the pairing energy that determines the magnitude of the phenomenological factorizable pairing 
potential of the form V(p, p') = — Vb<^(n)</>(n'), A is the amplitude of the superconducting order parameter 
A(p) = A(f>(n), n = p/p is a unit vector along the momentum, and the function </>(n) specifies the symmetry 
and anisotropy of A(p) in the momentum space [e. g., 0(n) = const for isotropic s-wave pairing; </>(n) = cos(2iy9) 
for a specific case of rf-wave pairing, where ip is an angle between the vector n and the 2>axis; </>(n) = r + cos(2<^) 
for a specific case of mixed (d+s)-wave pairing, the constant r being the measure of the partial weight of s-wave 
component in A(p)]. 

To find the dependence of V^ 1 on A, we make use of the set of mean-filed self-consistent equations for a 
superconductor containing both nonmagnetic and magnetic impurities (see Ref. [Jj): 
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where T is the temperature; £(p) is the quasiparticle energy measured from the chemical potential; u> = nT(2n+ 
1) are Matsubara frequencies; c n and c m are the concentrations of nonmagnetic and magnetic impurities, 
respectively; u n is the matrix element for potential electron scattering by an isolated nonmagnetic impurity; 

and v%H are the matrix elements for, respectively, potential and exchange (spin-flip) scattering by an 
isolated magnetic impurity. 

Restricting all electron momenta in Eqs. ©-© to the FS, replacing by N(0) J d£(k) J FS df2k/|<9£(k)/dk|, 
and integrating over £(k), one has 

V^ 1 =nN(0)(cj ) 2 (n)) FS Tj2f("^) > ( 6 ) 
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the angular brackets (...)fs stand for a FS average, 
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and we have introduced the electron relaxation times t p and r s for potential and spin-flip scattering, respectively: 



- = 2n(c n \u n \ 2 + c m |< ot | 2 )iV(0), - = 2nc m \u%\ 2 N(Q). 
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III. RESULTS 

Making use of a standard procedure, one can transform Eq. © into 
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where T c q is the value of T c in the absence of impurities. At A = 0, one obtains from Eqs. 0-©, ltl2ll an 
expression for T c as a function of potential and spin-flip scattering rates, p p = l/4irr p T c and p s = 1/2itt s T c , 
respectively (see Ref. jjj) 
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where p = p p + p s /2 is the total scattering rate and W is the digamma function. 

Expanding f(ui, A) in powers of A 2 up to A 2 and differentiating Eq. © with respect to A 2 , one has from 
Eq. © in the vicinity of T c 
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Next, taking Eq. 11 -til into account, one has from Eq. H12I1 an expression for A 2 (T) in the vicinity of T c 
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Here and below ip( n ' (z) are the polygamma functions [the n-th derivatives of the digamma function ip(z)] defined 

as V (n) ^) = (-l)" +1 "!Efco( fc + 

Substituting Eq. 11511 into Eq. 114L we obtain from Eq. the expression for the specific heat jump AC 

normalized by the specific heat in the normal state, 
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Finally, after simple but rather cumbersome calculations, we find from Eqs. 0-© an expression for the 
denominator in Eq. 111611 
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Eqs. ijlfijl and ifTTIl . together with Eq. l(T3|l for T c , give the value of the specific heat jump for a superconductor 
that is characterized by an arbitrary anisotropy of A(p) [i. e., by an arbitrary angular function <fi(n)] and 
contains, in general, both potential and spin-flip scatterers. In particular cases of (i) spin-flip scattering in 
an isotropic s-wave superconductor with <^(n) = const and (ii) potential scattering in a highly anisotropic 
(e. g., d-wave) superconductor with (0(n)) F s = 0, Eqs. 11 (it and II17|I reduce to the well-known expressions 
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respectively [22]. In the case of a superconductor that has an arbitrary anisotropy of A(p) but contains non- 
magnetic impurities only, Eqs. lfl6|l and ifTTIl reduce to the results of Harah et al. [22j . In the absence of any 
impurities (T c = T c q), one has 

AC 12 (d> 2 (n))l Q 



C„(T c0 ) 7C(3) <^(n)) FS ' 
1.202 is the Riemann zeta function. Eq. ll2fil) has been widely used to analyze the effect of 



where £(3) 

anisotropy of A(p) on the specific heat jump in clean superconductors (see, e. g., 23\). For a clean isotropic 
[4>(n) = const] superconductor we arrive at a familiar BCS result, AC / C n (T c0 ) = 12/7£(3) ~ 1.426. 



IV. DISCUSSION 

In what follows, we shall model the dependence of A(p) = A0(n) by the angular function cp(n) = r + 
cos(2(p). The value of r — corresponds to <i-wave pairing, while r — > 00 (Ar -^const) in an isotropic s-wave 
superconductor. The smaller is the value of r, the higher is the anisotropy of A(p). The moments (</>(n)) F s, 
(</> 2 (n)) FS , (0 3 (n)) FS , and (<?(> 4 (n)) F s that enter Eqs. O, and (H) are equal to r, r 2 + 1/2, r 3 + 3r/2, 

and r 4 + 3r 2 + 3/8, respectively. 

Note that the value of the specific heat jump in a clean superconductor is a nonmonotonous function of r. It 
follows from Eq. 12(H that the normalized specific heat jump initially decreases with r from AC / C n (T c o) w 0.951 
at r = down to AC / C n (T c o) ~ 0.666 at r = \/3/8 f» 0.612 and next increases again up to AC/C n (T c o) = 
12/7C(3) w 1.426 at r -> 00, see Fig. 1. 

Now let us analyze the behavior of AC /C n (T c ) versus T c upon addition of magnetic and/or nonmagnetic 
impurities to the initially clean sample with the critical temperature T c q. We note that nonmagnetic impurities 
result in the potential scattering only, while magnetic impurities generally result in both spin-flip and potential 
scattering. In this respect, the combined effect of nonmagnetic and magnetic impurities has much in common 
with the effect of magnetic impurities only, the difference being in the ratio of potential to spin-flip scattering 
rates as a function of impurity concentrations. 



A. Nonmagnetic disorder 



First, we consider the case that there are no magnetic impurities in the sample, i. e., p s = and, hence, 
p = p p (see also Ref. [0j). At low concentration of nonmagnetic impurities, i. e., at (T c q — T c )/T c q « 1, one 
has from Eqs. JOU and {UJ 
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where we took into account that p = (2/7r 2 )(l-T c /r cO )(0 2 (n)) F s/[(^ 2 (n)) FS -(0(n)) FS ] at (T c0 -T c )/T cQ « 1, 
see Eq. 1)1311 and Ref. [Jj. Note that the term in round brackets in Eq. 12111 changes sign from positive to negative 
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as r increases up to ro « 1.75. Analysis of Eqs. iflfijl . (tlTI) . and ll2l)l shows that in weakly anisotropic supercon- 
ductors with r > ro, the normalized specific heat jump increases monotonously up to AC / C n (T c ) = 12/7^(3) as 
T c is suppressed by nonmagnetic impurities, in a close agreement with the behaviour of AC / C n {T c ) in two-gap 
superconductors whose thermodynamics is, to some respect, similar to that of anisotropic superconductors, see 
Ref. [2I| (note that in an isotopic superconductor, r — > oo, the nonmagnetic disorder has no effect on both 
T c and AC/C n (T c ), see Refs. ^J|2£]). Contrary, at < r < r , the normalized specific heat jump initial- 
ly decreases with decreasing T c , passes through a minimum at T c /T c0 — (T c /T cQ )*, and then increases up to 
AC /C n {T c ) = 12/7£(3) as T c — > 0. The position of the minimum at the curve of AC / C n {T c ) versus T c /T c q 
depends on r, the value of (T c /T c q)* being reduced from 1 to as r decreases from ro down to zero [J]]. 

In a G?-wave superconductor without an admixture of s-wave, i. e., at r = 0, the normalized specific heat 
jump decreases monotonously down to zero as T c /T c q is suppressed from 1 to 0. For an arbitrary function </>(n) 
obeying the condition (</>(n))i?s = we have at T c /T c q « 1 



where 7 = e° w 1.781 and C is the Euler constant. Note that (</> 4 (n)) FS > {(f> 2 (n)) 2 FS for any 0(n), so the 
denominator in Eq. 112211 is always positive. For our choice of 4>(n) = r + cos(2</?) one has AC / C n {T c ) = 
(16/5h 2 (T c /T c0 ) 2 at r = and T c « T c0 . 

The effect of nonmagnetic impurities on the specific heat jump in superconductors with different anisotropy 
of the order parameter is illustrated in Fig. 2. The difference in the behavior of AC/C n (T c ) at T c — » in 
superconductors with r = and r ^ stems from the fact that, while T c of a d-wave superconductor vanishes 
at a finite value of p p — 1/47 ~ 0.140 (i. e., at a finite concentration of nonmagnetic impurities), the value of T c 
in a superconductor with the nonzero Fermi surface average of A(p) asymptotically goes to zero as p p increases 
at p s = 0, see Refs. [7ll26|. 



Although pure spin-flip scattering never happens in real materials since magnetic impurities give rise to not 
only spin-flip scattering but to potential scattering as well, we nevertheless (partly for pedagogical purposes) 
consider briefly the limiting case that a superconductor contains spin-flip scatterers only, i. e., p p = and 
p = p s /2. It is straightforward to show that in this case there are no qualitative differences among the curves 
of AC/C n (T c ) versus T c /T c q for different values of r. From Eqs. Illfit and l|17JI one finds that the normalized 
specific heat jump decreases monotonously down to zero with decreasing T c , no matter what the symmetry and 
the degree of anisotropy of A(p) are, see Fig. 3. The physical reason is that spin-flip scatterers, contrary to 
potential ones, are pair breakers in both d-wave and s-wave superconductors [Jj. Two quantitative differences 
between d-wave and s-wave symmetries are (i) the different values of AC/C n (T c ) at T c — T c0 and (ii) the 
different values of the coefficient f[(f>(n)] in the dependence AC/C n (T c ) — /[0(n)](T c /T c o) 2 at T c /T c q « 1. 
For a d-wave superconductor, the value of f[4>(n)] coincides with that in the case of nonmagnetic disorder, see 
Eq. i22\) . while / = 87 2 for an isotropic s-wave superconductor, in agreement with Ref. |12j . 



Now we turn to a general case that there are both nonmagnetic and magnetic impurities in a superconductor. 
To begin with, we consider a situation when purely nonmagnetic impurities are added to a superconductor 
that already contains a small quantity of magnetic impurities and, as a consequence, initially has the critical 
temperature T^ lower than the value of T c $ in the absence of any impurities. The higher is the concentration 
of magnetic impurities, c m , the greater is the value of 5T co /T co , where ST co — T c o — T' c0 . We assume that the 
value of c m remains unchanged upon increase in the concentration of nonmagnetic impurities and corresponding 
decrease of T c , i. e., p s o = 1/2ttt s T c o = const. Since our prime interest here is with the case ST co /T co « 1, we 
plot AC / C n (T c ) versus T c /T c q rather than T c /T c ' . 

In a d-w&ve superconductor, i. e., at r = 0, the dependence of AC/C n (T c ) on T c /T c q for any value of 
ST co /T co is the same as in the absence of spin-flip scattering, see Fig. 2. This is because at {(/)(n)) F g = both 
T c and AC are functions of the total scattering rate p = p p + p s /2 only, irrespective of the scatterers' type, 
see Eqs. l(T3ll . ltl6)l . and ifTTIl . Fig. 4 shows the curves of AC /C n (T c ) versus T c /T c0 for different values of r 7^ 
and ST c q/T c q. At low r = 0.2, i. e., in a [d + s)-wave superconductor with a small admixture of s-wave, or 
in a strongly anisotropic s-wave superconductor, the presence of finite, though very small (ST c q/T c q < 0.01) 
amount of spin-flip scatterers in the sample results in a drastic change in the dependence of AC /C n (T c ) on 




(22) 



B. Pure spin-flip disorder 



C. Combined nonmagnetic and magnetic disorder 

1. Constant concentration of spin-flip scatterers 
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T c /T c o. After passing through a minimum, AC/C n (T c ) does not increase up to 12/7£(3), as in the absence of 
magnetic impurities, but reaches the maximum and then decreases again down to zero as T c — > 0, see Fig. 4. It 
is straightforward to show from Eqs. i|13ll . (IKill . and II17|I that 

AC 1 (T c \ 2 _tt 4 
C n (T c ) ~ 2p 2 s0 \T c0 ) "32 

at T c << ST c q. As the concentration of spin-flip scatterers (and, hence, the value of ST c q/T c q) increases, the 
maximum of AC/C n (T c ) decreases in height and gradually disappears. Starting with ST c0 /T cQ w 0.02 there are 
neither minimum nor maximum of AC/C n (T c ) versus T c /T c0 curve. We note that for T c0 « 100 K such values 
of 8T c0 /T c0 correspond to commonly observed small (several Kelvins) variations of T c0 between the samples 
obtained under slightly different conditions that can reflect the different amount of magnetic impurities in the 
samples. 

As r increases, i. e., as the anisotropy of A(p) becomes weaker, the tendency in the change of the specific 
heat jump upon increase in the initial concentration of magnetic impurities is qualitatively conserved, see Fig. 
4. Quantitatively, an increase in r results in the increase in the value of ST c0 /T c o above which AC/C n (T c ) 
becomes a featureless monotonous function of T c /T c q. In particular, both the minimum and the maximum of 
AC / C n (T c ) versus T c /T c q curve disappear at 5T c q/T c q ps 0.1, 0.25, and 0.4 for r = 0.6, 1, and 2, respectively. So, 
the sensitivity of the specific heat jump to magnetic impurities is higher in strongly anisotropic superconductors 
with small but nonzero values of r. 



(0 2 (n))F5 
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2. Constant ratio of spin-flip to potential scattering rates 

Let us now consider the case that the relative contribution from spin-flip scattering to the total scattering 
rate, p = p p + p s /2, remains constant upon disordering, i. e., the value of the dimensionless coefficient 

a - T " - (24 ) 

Tp X +r7 l (c„/c m )| U „| 2 + |^ ot | 2 + |<-| 2 Pp + Ps/2' 

see Eq. l(TT|l . does not change upon addition of magnetic (and, in general, nonmagnetic) impurities. This holds, 
first, if a superconductor is doped by magnetic impurities only (i. e., c„ = 0) and, second, if the ratio of 
nonmagnetic to magnetic impurity concentrations, c„/c m , remains unchanged, see Eq. Il24t . The latter is a 
reasonable approximation for doping by given chemical elements or irradiation by a given type of particles, at 
least at relatively low (but sufficient to destroy the superconductivity) doping levels or radiation doses. 

Thus for a given degree of A(p) anisotropy (i. e., in our model, for a given value of r), the dependence of 
AC / C n (T c ) on T c /T c o is governed by the value of material-dependent and disorder-dependent coefficient a. 
The greater is the relative contribution from exchange scattering by magnetic impurities to the total scattering 
rate, the higher is the value of a. In general, a ranges from in the absence of exchange scattering to 1 in the 
absence of potential scattering. Note, however, that since there always exist two channels of carrier scattering 
by magnetic impurities (potential and spin-flip ones), see Eq. lUTf . the value of a is less than unity even at 
c„ = 0. Below we consider the case a « 1 that seems to be relevant to the experimental situation. 

It follows from Eqs. l(T3ll . ijlfil) . and ifTTj) that in a d-wave superconductor with r = the curves of AC /C n {T c ) 
versus T c /T c q are the same for any value of a in the whole range of a, see Fig. 2 and Fig. 3. Contrary, at r ^ 
the specific heat jump appears to be extremely sensitive to spin-flip scattering of charge carriers. Fig. 5a shows 
the dependencies of AC / C n (T c ) on T c /T c q in a strongly anisotropic non d-wave superconductor with r = 0.2 
for different values of a. One can see that increase in a results in a gradual disappearance of the minimum 
(and maximum) of AC/C n (T c ). In the presence of even a minor spin-flip component in the scattering potential, 
a w 0.02, the normalized specific heat jump decreases monotonously as T c decreases from T c q down to zero, 
and the curve of AC / C n (T c ) versus T c /T c q looks like that in a d-wave superconductor. As the gap anisotropy 
weakens (i. e., the value of r increases) the "critical "value of a above which the normalized specific heat jump 
starts to decrease monotonously under disordering first increases up to ss 0.05 atr» 1 and next decreases 
again, see Fig. 5. 



D. Implications for the experiment 

Numerous experiments on various superconductors, including borocarbides Yi_ a; R a ;Ni2B2C (R = Gd, Dy, 
Ho, Er) pl|, organic compound (TMTSF) 2 C10 4 [23, Ui-^Th^Bea [2E], HTSCs YBa2(Cu 1 _ x M x ) 3 07-<5 (M = 
Zn [HI EDJ EH E2, Fe ESI' Ni EB' Cr E2) and Lai.85Sro.i5Cui_ v Zn J/ 04 ESI, etc have revealed that the value 
of AC / C n (T c ) decreases monotonously as T c is suppressed by impurities. To the best of our knowledge, there 
were no experimental indications for the nonmonotonous behaviour of AC / C n (T c ) in disordered superconduc- 
tors. Note, however, that the chemical substitution results not only in the suppression of T c and decrease of 
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AC/C„(T C ) but also in a very strong broadening of the superconducting transition. As a consequence, the spe- 
cific heat anomaly is rapidly smeared out by the disorder, so that the dependence of AC '/ C n (T c ) on T c /T c q can 
be determined more or less reliably, in the best case, at T c /T c q > 0.3 ~ 0.4 only. Meanwhile, it follows from the 
results presented above that the value of T c /T c q below which AC/C n (T c ) starts to increase under disordering 
depends on the degree of A(p) anisotropy and is very small in strongly anisotropic superconductors. 

Recently Zhao has fitted A(p) to single-particle tunneling and angle-resolved photoemission spectra of 
YBa2Cu307-«5 ; 2]. To compare his fit with our model form of A(p), it is convenient to introduce the coefficient 
X = 1 — (A(p))|, s /(A 2 (p)}i?s as a measure of the degree of in-plane anisotropy of A(p), where (...)fs means 
the Fermi surface average. The range < \ < 1 covers the cases of isotropic s-wave (A(p)=const, \ = 0), 
d-wave ((A(p))_fs = 0, x = 1); an d mixed (d + s)-wave or anisotropic s-wave (0 < x < 1) symmetries of A(p). 
Making use of the results presented in Ref.p], one has x ~ 0.9 for YBa2Cu3C>7-<5- Since x = 1/(1 + 2r 2 ) for 
A(p) = A[r + cos(2(p)], this value of \ corresponds to r w 0.2. We note that the choice of x ~ 0.9 allows 
for a quantitative explanation 34J of the quasilinear decrease of T c in electron-irradiated YBa2Cu307_5 single 
crystals As follows from Figs. 4 and 5, at r = 0.2, the upturn of AC / C n (T c ) takes place at T c /T c q w 0.1, 
either in the absence or at a very low concentration of spin-flip scatterers. 

So, the necessary condition for the AC /C n {T c ) upturn at low T c /T c0 is, except for the non-pure d-wave 
symmetry of A(p), a relatively small contribution of spin-flip scattering to the total scattering rate. Besides, a 
superconductor should be disordered very uniformly so that the transition width AT C remained lower than T c 
down to as low as possible T c values, in order to preserve a clear specific heat anomaly at T c and to make possible 
the experimental determination of the AC / 'C n (T c ) versus T c /T c q curve in a wide region of T c /T c q values. In 
this respect, the irradiation-induced disorder has advantages over the chemical substitution. For example, in a 
recent paper |82] , Rullier- Albenque et al. reported the results of experimental studies of T c degradation under 
electron irradiation of underdoped and optimally doped YBa2Cu307_5 single crystals. The authors of Ref. jHI] 
succeeded in creation of an extremely uniform distribution of radiation defects over the sample, so that the 
value of AT C never exceeded 5 K. Moreover, the value of AT C did not increase monotonously with radiation 
dose but had a maximum at T c /T c0 « 0.3 and next decreased down to AT C < 1 K at the highest dose for which 
the resistive superconducting transition still was observed at T c w 1 K. According to the theoretical fit [331 1° 
the experimental data |35,, at x = 0.9, the value of a = 0.01 ± 0.01 is low enough for the upturn of AC/C n (T c ) 
be observable at T c /T c q « 0.1, see Fig. 5. Hence, it is of great interest to study the behavior of AC versus T c /T c q 
in such samples. 

Finally, we note that the experimentally observed nonuniversality of AC / C n (T c ) versus T c /T c q curve that 
has been previously ascribed to the carrier concentration effects [2lj may in fact be (at least partly) due to 
different contributions of spin-flip scattering to pair breaking in different superconducting materials and / or for 
different doping elements. We note also that it would be very interesting to study experimentally the behavior 
of AC/C n (T c ) versus T c /T c o in non-cuprate superconductors with different degree of the superconducting gap 
anisotropy and various ratios of spin-flip to potential scattering rates. 

V. CONCLUSIONS 

We have shown that in a pure d-wave superconductor, the normalized specific heat jump AC/C n {T c ) 
decreases monotonously upon disordering by both nonmagnetic and magnetic defects or impurities. So, in d- 
wave superconductors, AC /C„(T C ) is, to a first approximation (keeping in mind the assumptions made), a 
universal function of T c /T c0 , i. e., it does not depend on the relative contribution of spin-flip scattering to 
the total scattering rate and, hence, on a specific type of defects and impurities. On the other hand, under 
nonmagnetic disordering of a superconductor with a nonzero Fermi surface average of the order parameter, 
AC /C n (T c ) initially decreases with decreasing T c , passes through a minimum and then increases again. The 
minimum at the curve of AC/C n (T c ) versus T c /T c q moves to higher values of T c /T c o as the anisotropy of the 
order parameter becomes weaker. 

In disordered strongly anisotropic non ci-wave superconductors, AC /C n {T c ) is extremely sensitive to spin- 
flip scattering of charge carriers. At relatively weak spin-flip scattering, AC/C n (T c ) becomes a featureless 
monotonous function of T c /T c q. So, the spin-flip scattering of charge carriers removes the qualitative difference 
between the dependencies of AC '/C n (T c ) versus T c /T c0 in disordered d-wave and (d + s)-wave (or anisotropic 
s-wave) superconductors. Hence, it would be very difficult to discriminate between pure <i-wave and non pure 
d-wave A(p) if the concentration of spin-flip scatterers in the sample is higher than a certain critical value. So, 
to observe the nonmonotonous dependence of AC / C n (T c ) on T c /T c q in anisotropic superconductors, one should 
(i) make use of uniformly disordered samples with a clearly pronounced specific heat anomaly at T c down to 
low T c /T c q values and (ii) minimize the concentration of the spin-flip scatterers in the sample. 
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Fig. 1. Specific heat jump AC normalized by the normal state specific heat C„(T c o) versus the coefficient 
r that specifies the anisotropy of the superconducting order parameter A(p) = A[r + cos(2(^)], for a clean 
superconductor without any impurities. Dashed line shows the value of AC/C n (T c o) = 12/7^(3) ~ 1.426 in an 
isotropic s-wave superconductor. 
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T c lT c o 

Fig. 2. Specific heat jump AC normalized by the normal state specific heat C n (T c ) versus the normalized 
critical temperature T c /T c q for a superconductor disordered by nonmagnetic impurities. The superconducting 
order parameter is assumed to have the form A(p) = A[r + cos(2<p)], where r = (closed circles); 0.2 (squares); 
0.6 (triangles); 1 (open circles); 2 (pluses). Solid lines are guides for the eye. Dashed line corresponds to an 
isotropic s-wave (r — > oo, Ar — >-const). 
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Fig. 4a. Specific heat jump AC normalized by the normal state specific heat C„(T C ) versus the normalized 
critical temperature T c /T c q for a superconductor with the order parameter A(p) = A[r + cos(2</?)]. The super- 
conductor initially contains a small amount of spin-flip scatterers and thus has the initial critical temperature 
T' cQ < T c q. It is further disordered by nonmagnetic impurities only, so that the concentration of potential scat- 
terers exceeds that of spin-flip scatterers. r — 0.2. ST co /T c q = (T c q — T' c0 )/T c $ = (solid line), 0.001 (pluses), 
0.003 (triangles), 0.01 (squares), and 0.03 (circles). These values of 5T co /T c q correspond to the values of the 
spin-flip pair breaking rate p s0 = 1/2ttt s T c0 = 0, 0.00038, 0.00113, 0.00376, and 0.01127 respectively. Note 
that at T c /T c0 > 0.2 the curves of AC / C n {T c ) versus T c /T c0 for different p s o (i. e., for different 5T co ) almost 
coincide. 
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Fig. 4b. Specific heat jump AC normalized by the normal state specific heat C n (T c ) versus the normalized 
critical temperature T c /T c q for a superconductor with the order parameter A(p) = A[r + cos(2<p)]. The super- 
conductor initially contains a small amount of spin-flip scatterers and thus has the initial critical temperature 
T^Q < T c q. It is further disordered by nonmagnetic impurities only, so that the concentration of potential scat- 
terers exceeds that of spin-flip scatterers. r = 0.6. 5T co /T c q = (solid line), 0.01 (pluses), 0.03 (triangles), 0.1 
(squares), and 0.3 (circles). [p s0 = 0, 0.00285, 0.00854, 0.02823, and 0.08232 respectively]. 
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Fig. 4c. Specific heat jump AC normalized by the normal state specific heat C n (T c ) versus the normalized 
critical temperature T c /T c o for a superconductor with the order parameter A(p) — A[r + cos(2</?)]. The super- 
conductor initially contains a small amount of spin-flip scatterers and thus has the initial critical temperature 
T^Q < T c q. It is further disordered by nonmagnetic impurities only, so that the concentration of potential scat- 
terers exceeds that of spin-flip scatterers. r = 1. 5T co /T c q = (solid line), 0.03 (pluses), 0.1 (triangles), 0.2 
(squares), and 0.3 (circles). [p s0 = 0, 0.00727, 0.02400, 0.04728, and 0.06975 respectively]. 
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Fig. 4d. Specific heat jump AC normalized by the normal state specific heat C n (T c ) versus the normalized 
critical temperature T c /T c q for a superconductor with the order parameter A(p) = A[r + cos(2<p)]. The super- 
conductor initially contains a small amount of spin-flip scatterers and thus has the initial critical temperature 
T^Q < T c q. It is further disordered by nonmagnetic impurities only, so that the concentration of potential scat- 
terers exceeds that of spin-flip scatterers. (d) r = 2. ST CO /T c o = (solid line), 0.03 (pluses), 0.1 (triangles), 0.3 
(squares), and 0.5 (circles). [p s0 = 0, 0.00641, 0.02114, 0.06116, and 0.09731 respectively]. 
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Fig. 5a. Same as in Fig. 4 for a superconductor with a constant ratio a of spin-flip to potential scattering 
rates, r = 0.2. a = (solid line), 0.001 (pluses), 0.003 (triangles), 0.01 (squares), and 0.03 (circles). Note that 
at T c /T c o > 0.2 the curves of AC/C n (T c ) versus T c /T c q for different a almost coincide. 
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Fig. 5d. Same as in Fig. 4 for a superconductor with a constant ratio a of spin-flip to potential scattering 
rates, r = 2. a = (solid line), 0.001 (pluses), 0.01 (triangles), 0.03 (squares), and 0.1 (circles). 
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